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The modern theory of neutrinoless double beta decay includes a scaling factor that has often 
been treated inconsistently in the literature. The nuclear contribution to the decay half life can be 
suppressed by 15-20% when scaling factors are mismatched. Correspondingly, (rui/) is overestimated. 

PACS numbers: 


I. INTRODUCTION 

In recent years, experimental evidence for neutrino masses and mixing have lead to a concerted effort to refine the 
methods used in calculating reaction rates for double-beta decay. One of the largest uncertainties in these calculations 
is the determination of the nuclear matrix elements (NME). Together with kinematic factors and experimental bounds 
on the decay half life, an average neutrino mass can be extracted. An accurate determination of the NME is crucial 
and improvements in the QRPA and shell model techniques used to calculate them continue to be explored. 

In the modern theory of neutrinoless double beta decay (/3/3(Oi^)) definitions of the NME include a scaling factor 
introduced such that the NME is dimensionless. When using these scaled NME to determine the /3/3(0^) decay 
rate, the scaling factor must be compensated for elsewhere. However, the scaling factor has not always been treated 
consistently in the literature. In this article, we detail how the scaling factor is introduced into the theory and 
illustrate the 15-20% suppression of the nuclear contribution when mismatched scaling factors are used. 


II. NEUTRINOLESS DOUBLE BETA DECAY 


In the simplest form of the weak Hamiltonian 


Hw 



the half life of the 0+ ^ 0'^' l3(3{Qv) is written as 
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where is the NME, Gqi is the so-called phase space or kinematic factor and nii, and We are the neutrino and 

electron masses respectively. The NME is given by 

MfO") = (3) 
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The sum is taken over intermediate states Na- A detailed explanation of the NME can be found in P|. We will use 
notation adapted from Q throughout this article. 

The v potential induced by the virtual v exchange is given by 
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The initial, intermediate and final state energies are denoted by Ei, Ea and Ef respectively and Ci denotes the energy 
of the ith electron. 

The kinematic factor Gqi in Eq. © is defined as 
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We have assumed S' = 0 electron wave functions with no r dependence. The Fermi functions, Eo(2', e), depend upon 
the charge of the daughter nucleus, Z, and the energy of the zth electron, a. 

In the early (3l3{0v) calculations, the NME and kinematic factor were defined as above. The NME were given in 
units of fm“^ and the kinematic factors in units of yr“^ fm^ (for example, 0 )- In the mid-eighties, a scaling factor 
was introduced into the /3/3(Oi/) theory. The v potential, Eq. (TJ, is scaled by a factor of i? = such that the 

NME are dimensionless: 
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The scaling factor, ii, in h-^. is compensated for by introducing 1/R^ into the definition of the kinematic factor, Eq. 

0: 

G^oiiR) = 1^2 / 

The Gqi{R) have been calculated by many authors; for example, 0 using R = 1.2A^/^ and 0 using R = lAA^/^. 
Though the underlying physics of the kinematic factor is unchanged, the published values of Gqi (R) are significantly 
different due to different choices of R. Provided that the NME and kinematic factor are calculated using the same 
scaling factor, these differences are irrelevant. However, the R used in calculating the NME have not always been 
carried consistently to the kinematic factor. 

In recent years, numerous calculations of NME have been performed using h^. As techniques develop, it is customary 
to compare with previously published values. These comparisons are often made by defining Gqo: 
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If scaling is treated consistently. Coo is independent of the scaling factor. However, in citing previous calculations of 
the NME, the scaling factor has not always been accounted for properly. Often, Gqi{R = 1 . 2 A^C) jg used to calculate 
Coo regardless of the scaling factor used in determining the NME. If one combines NME calculated with R = 1.1 
fm with Gqi{R = 1 . 2 A^C)^ is suppressed by (1.1/1.2)^ ~ 20%. Correspondingly, this leads to an overestimation 
of {rriy). 

Tablejp includes several Coo^predictions for ^®Ge /3/3(Oi^). This table was originally published in a recent review 
article 0 and adapted from The mismatch of scaling factors has occurred several times in the literature. We 
include revisions to Table 2 of 0 because it is one of the more thorough listings of NME calculations. 

In Table U we include both the previously published Gqo and the revised values obtained using the correctly scaled 
GQi(i?). For clarity, the value of tq used in the original publication is included. In some instances rp was not clearly 
stated; the rg value was extracted from T 1 / 2 , (m^) or stated Gqi{R) values given in the original paper assuming that 
the scaling was originally treated consistently. 






3 


Using the simple S = 0 electron wave functions with no r dependence, we obtain Gqi{R = 1.2A^/^) = 6.46 • 10 
yr-i and G^{R = = 7.78 • IQ-^^ yr ^. These values are within a few percent of those published previously 

by 0 and j3|. Using the appropriate kinematic factor, most Coo determined using NME with tq = 1.1 are changed 
by « 20%. Assuming a half life of 4 • lO^"^ yr the spread of predicted {rriy) values is unchanged, 0.022-0.068. However, 
several predicted neutrino masses are reduced. 


/^revised 

^00 

X 10 - i ^ yr - i ) 

/^old 

^00 

X 10 - i ^ yr - i ) 

(m^) 

(eV) 

Method 

ro 

(fm) 

Reference 

11.2 

11.2 

0.024 

QRPA 

1 . 2 ^ 


6.97 

6.97 

0.031 

QRPA 

NS 


7.51 

7.51 

0.029 

Number-projected QRPA 

NS 


7.19 

7.33 

0.030 

QRPA 

1.1 

[ 3 j 

12.1 

12.0 

0.023 

QRPA 

NS 


13.3 

13.3 

0.022 

QRPA 

NS 

[Ifl] 

8.34 

8.27 

0.028 

QRPA 

1.2 

liJJ 

1.89 - 12.8 

1 . 85 - 12.5 

0 . 059 - 0.023 

QRPA 

1 . 2 ^ 

[ 12 j 

5.02 - 5.93 

1 . 8 - 2.2 

0 . 036 - 0.033 

QRPA 

l.U 

lia 

8.61 

8.36 

0.028 

QRPA 

1 . 2 ^ 


1.40 

1.42 

0.068 

QRPA with up pairing 

1.1 


5.59 

4.53 

0.034 

QRPA with forbidden 

l.U 

M 

10.1 

10.3 

0.025 

RQRPA 

1.1 


6.10 

6.19 

0.033 

RQRPA with forbidden 

1.1 

ua 

6.77 - 7.72 

5 . 5 - 6.3 

0 . 031 - 0.029 

RQRPA 

i.U 

]m 

2.26 - 9.04 

2 . 21 - 8.83 

0 . 054 - 0.027 

RQRPA 

1 . 2 ^ 

[ 12 j 

4.48 

3.63 

0.038 

RQRPA with forbidden 

i.U 

M 

2.71 

2.75 

0.049 

Full RQRPA 

1.1 

m 

3.72 - 8.75 

3 . 36 - 8.54 

0 . 042 - 0.027 

Full RQRPA 

1 . 2 ^ 

[ 12 j 

6.66 - 9.43 

6 . 50 - 9.21 

0 . 031 - 0.026 

Second QRPA 

1 . 2 ^ 

[ 12 j 

0.34 - 0.40 

0 . 27 - 0.32 

0 . 139 - 0.128 

Self-consistent QRPA* 

l.U 

]m 

28.8 

28.8 

0.015 

VAMPIR* 

NS 

m 

15.6 

15.8 

0.020 

Shell model truncation* 

NS 

m 

7.03 - 16.2 

6 . 87 - 15.7 

0 . 030 - 0.020 

Shell model truncation* 

1 . 2 ^ 

[ 18 ] 

1.94 

1.90 

0.058 

Large-scale shell model 

1 . 2 ^ 

[ 19 ] 


TABLE I: A comparison of the nuclear contributions to T 1/2 , Coo, and the resulting {m^) values assuming T \/2 = 4.0 • 10^^ 
yr. Included are the values of ro used in each NME calculation. The ro has been explicitly stated in the original publication 
unless otherwise indicated as: NS (no scaling/Coo published in original paper); ^ (ro inferred from published values of (rn^), 
ri /2 or Goi) and ^ (Published Coo; ro inferred from published values of {m^), T 1/2 or Coi) 

It is important to point out Gqi, defined without the scaling factor, does depend on R through the Fermi functions 
of the electron wave functions. The appropriate choice of R depends upon the nucleus being considered and should be 
chosen such that experimental values of the mean square radius (r^) are reproduced. For a uniform charge distribution, 
r2 = 5/3(r2). In TableUD we give the unsealed Gqi calculated using R values fit to experimental root-mean-square 
radii when possible Comparing Gqi for it is clear that the unsealed Gqi are not very sensitive to the 

choice of R. The significant differences obtained by 0 and 0 for Gq^ (i?) are predominately due to the scaling factor. 

To avoid confusion in the future, we strongly encourage that further calculations of NME be published with either 
no scaling factor included, or a clear indication of what choice the authors have made for tq. 
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Nucleus 

(^2)1/2 

fm 

ro 

fm 

/-lOv 

^01 

fm^ yr-** 

00 

Q 

3.470 

1.23 

1.2362;10"***^ 


3.451 

1.23 

1.2362;10"***^ 

"®Ge 

4.081 

1.24 

1.6632;10"***^ 

»2Se 

no data 

1.20 

7.784x10“***^ 

^®Zr 

4.396 

1.24 

1.792a:10“***^ 

^°°Mo 

4.430 

1.23 

1.436a:10“***^ 

“®Cd 

4.639 

1.23 

1.7202:10“***^ 

128 Te 

no data 

1.20 

e.oooxio"**-* 

isoTe 

no data 

1.20 

1.6612:10“***^ 

i36Xe 

no data 

1.20 

1.8252:10“***^ 

isoNd 

5.048 

1.23 

8.7192:10"***^ 


5.015 

1.22 

8.7502;10"***^ 


4.948 

1.20 

8.8132:10"***^ 


TABLE II: 0*^ —> O’*" PP{Ou) kinematic factor Goi calculated at the specified values of ro obtained by fitting the experimental 
when available 1 ^ . 
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